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Abstract. Context. The nature of the progenitors of type la supemovae is still under controversial debate. KPD 1930+2752 is 
one of the best SN la progenitor candidates known today. The object is a double degenerate system consisting of a subluminous 
B star (sdB) and a massive white dwarf (WD). Maxted et al. 1 2000^1 conclude that the system mass exceeds the Chandrasekhar 
mass. This conclusion, however, rests on the assumption that the sdB mass is 0.5 Mq. However, recent binary population syn- 
thesis calculations suggest that the mass of an sdB star may range from 0.3 Mq to more than 0.7 Mq. 

Aims. It is therefore important to measure the mass of the sdB star simultaneously with that of the white dwarf. Since the 
rotation of the sdB star is tidally locked to the orbit the inclination of the system can be constrained if the sdB radius and the 
projected rotational velocity can be measured with high precision. An analysis of the ellipsoidal variations in the light curve 
allows to tighten the constraints derived from spectroscopy. 

Methods. We derive the mass-radius relation for the sdB star from a quantitative spectral analysis of 150 low-resolution spectra 
obtained with the Calar Alto 2.2 m telescope using metal-rich, line blanketed LTE model atmospheres with and without NLTE 
line formation. The projected rotational velocity is determined for the first time from 200 high-resolution spectra obtained with 
the Keck I 10 m, and the ESO-VLT 8.2 m telescopes. In addition a reanalysis of the published light curve is performed. 
Results. The atmospheric and orbital parameters are measured with unprecedented accuracy. In particular the projected rota- 
tional velocity Vi„, sin / = 92.3 + 1 .5 km s" ' is determined. Assuming the companion to be a white dwarf, the mass of the sdB is 
limited between 0.45 Mq and 0.64 Mq and the corresponding total mass of the system ranges from 1.33 Mq to 2.04 Mq. This 
constrains the inclination to / > 68°. The photometric analysis allows to constrain the parameters even more. A neutron star 
companion can be ruled out and the mass of the sdB is limited between 0.45 Mq and 0.52 Mq. The total mass of the system 
ranges from 1.36 Mq to 1.48 Mq and hence is likely to exceed the Chandrasekhar mass. The inclination angle is 80° and the 
light curve shows weak and shallow signs of eclipses. A high precision light curve is needed in order to accurately measure 
these eclipses. So KPD 1930+2752 qualifies as an excellent double degenerate supernova la progenitor candidate. 
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general 
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1. Introduction 

Supernovae of type la (SN la) play an important role in the 
study of cosmic evolution. They are regarded as the best stan- 
dard candles for the determination of the cosmological param- 
eters Ho, O. and A (e.g. Riess et al. .1998. Leibundgut 200L 
Perlmutter et al. I1999> together with the WMAP probe for 
anisotropics in the cosmic background radiation (e. g. Bennett 
et al. 2003 ). Although some of the most important discoveries 
in modern cosmology, in particular the accelerated expansion 
of the universe, were initially derived from distance measure- 
ments of SN la with high redshifts, the nature of their progen- 
itors is still under debate (Livio 2000 1. The progenitor popu- 
lation is a crucial information required to back the assumption 
that distant SN la can be used as standard candles as it was 
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done with the ones in the local universe. There is general con- 
sensus that only the thermonuclear explosion of a white dwarf 
(WD) is compatible with the observed features of SN la. For 
this a white dwarf has to accrete mass from a close compan- 
ion to reach the Chandrasekhar limit of 1.4 Mq (Hamada & 
Salpeter'196ri. Two main scenarios of mass transfer are cur- 
rently under discussion. In the so-called single degenerate (SD) 
scenario (Whelan & Iben ll973t . the mass donating component 
is a red giant/subgiant, which fills its Roche lobe and is contin- 
ually transfering mass to the white dwarf. According to the so- 
called double degenerate (DD) scenario (Iben & Tutukov 1984) 
the mass donating companion is a white dwarf, which eventu- 
ally merges with the primary due to orbital shrinkage caused 
by gravitational wave radiation. 

A progenitor candidate for the DD scenario must have a 
total mass near or above the Chandrasekhar limit and has to 
merge in less than a Hubble time. Systematic radial velocity 
(RV) searches for DDs have been undertaken (e.g. Napiwotzki 
l2U03 'and references therein). The largest of these projects was 
the ESO SN la Progenitor Survey (SPY). More than 1 000 WDs 
were checked for RV-variations (e. g. Napiwotzki et al.^OOS I. 
SPY detected ~ 100 new DDs (only 18 were known before). 
One of them may fulfill the criteria for progenitor candidates, 
even though the error margins are quite high (Napiwotzki et al. 
120051 Karl lMl^l . 

KPD 1930-H2752 was identified as a subdwarf B (sdB) 
star in the Kitt Peak - Downes survey (Downes 1986 1. The 
parameters which were derived from spectroscopy by model 
atmosphere fits (Saffer & Liebert J 995 J are consistent with 
the theoretical instability strip predicted by Charpinet et al. 
J1996> . After the first pulsating sdB stars with short periods 
(EC 14026 stars) were discovered (Kilkenny et al. 1997; Koen 
et al. nWTI Stobie et al HWl . Billeres et al. (2000) initiated 
a survey to search for such objects in the northern hemisphere. 
KPD 1930+2752 was selected from the list of Saffer & Liebert 
J1995> for their fast photometry program. They detected mul- 
tiperiodic variations with short periods and low amplitudes. In 
addition to 44 p-mode pulsations they found a strong variation 
at a much longer period of about 4100 s. This variation could 
be identified as an ellipsoidal deformation of the sdB star most 
likely caused by a massive companion. Billeres et al. (2000) 
predicted the period of the binary to be two times the period of 
this brightness variation (P = 8217.8 s = 0.0951 1 1 d). 

This was proven by Maxted et al. J2000> . who measured a 
radial velocity curve of KPD 1930+2752 which matched the 
proper period. The radial velocity amplitude K - 349.3 + 
2.7 km s"' combined with an assumption of the canonical mass 
for sdB stars MsdB = 0.5 led to a lower limit for the mass of 
the system derived from the mass function. This lower limit 
M > 1 .47 Mq exceeded the Chandrasekhar mass of 1 .4 Mq 
(Hamada & Salpeter 1961 1. Because there was no sign of a 
companion in the spectra, it was concluded that the unseen ob- 
ject must be a white dwarf. The system is considered to become 
double degenerate when the subdwarf eventually evolves to a 
white dwarf. Orbital shrinkage caused by gravitational wave 
radiation will lead to a merger of the binary in about 200 Myr 
which is much less than a Hubble time. Combining all the evi- 



dence, Maxted et al. concluded that KPD 1930+2752 could be 
a good candidate for the progenitor of a Type la supernova. 

From the theoretical point of view Ergma et al. J2001> 
questioned the double degenerate scenario in the case of 
KPD 1930+2752. Their simulations based on the derived pa- 
rameters of Maxted et al. ( 2000 1 suggested that the merger will 
lead to a single massive ONeMg white dwarf rather than a su- 
pernova explosion. The merger of two white dwarfs is a rather 
complicated process compared to slow accretion of material 
onto a heavy white dwarf. Detailed merger simulations in 3D 
are not yet avaible, but urgently needed to learn more about this 
possible SN la progenitor channel. 

The subluminous B stars form the hot end of the horizontal 
branch, the so called Extreme Horizontal Branch (EHB, Heber 
11986). The EHB models imply that they are core helium burn- 
ing objects with tiny hydrogen dominated envelopes, too thin 
to sustain hydrogen burning. The core mass is fixed by the on- 
set of the core helium flash at the tip of the red giant branch 
and depends only slightly on metalhcity and helium abundance. 
Hence the canonical core mass is resticted to a very narrow 
range of 0.46 to 0.5 Mq. Unlike normal horizontal branch stars, 
EHB stars do not ascend the asymptotic giant branch (AGB) 
after core helium burning has terminated, but evolve to higher 
temperatures until the white dwarf cooling sequence is reached. 

Binary evolution scenarios have been proposed (e.g. Han 
et al., 20d2| I2003> to explain the large fraction of sdB stars 
in binaries. These scenarios involve strong mass-transfer either 
through a common-envelope ejection, or through stable Roche 
lobe overflow. A merger of two helium core white dwarfs is 
another viable option. 

The binary population synthesis models of Han et al. <2002> 
predict a mass distribution for sdB stars that peaks near the 
canonical mass, but is much wider, ranging from 0.3 Mq to 
more than 0.7 Mq. The lower limit is based on the fact that no 
core helium burning can be sustained for stars with masses less 
than 0.3 Mq. The highest masses result from mergers. Hence 
it might be premature to assume that the mass of an sdB star 
is half a solar mass, as has been done in previous studies. 
Therefore we drop this assumption and treat the mass of the 
sdB component as a free parameter. 

Another drawback of all previous investigations was the 
lack of information on the inclination angle. We aim at deriv- 
ing constraints on the inclination for the first time by means 
of accurate measurements of the projected rotational velocity 
and the surface gravity from appropriate high-quality spectra 
using sophisticated model atmospheres. Because the rotation 
of the sdB star is tidally locked to the orbit, we can derive sin / 
as a function of the sdB mass from these two quantities. This 
approach combined with a reanalysis of the published light 
curve allows to constrain the masses of both components and, 
therefore, to check whether or not KPD 1930+2752 qualifies 
as a SN la progenitor candidate through a double degenerate 
merger. 

Section|2describes the newly obtained spectra, from which 
an improved radial velocity curve has been derived (Sect. |3}- 
For the quantitative spectral analysis, existing grids of metal 
line blanketed LTE models have been used and a new (hy- 
brid) approach is applied that allows to treat both departures 
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Fig. 1. Radial velocity curve. The data points are from the new 
HIRES spectra, the solid line is the fitted sine curve to these 
points. The residuals in the lower part show no signs of eccen- 
tricity. 

from LTE as well as metal line blanketing (Sect.|3i. Combining 
these results, described in Sect.|5l with the projected rotational 
velocity (Sect. |5) we constrain masses and inclination of the 
KPD 1930+2752 system (Sect.Q. Within the derived pai-ame- 
ter space synthetic light curves were modelled and fitted to the 
published light curve (Sect.|S} in order to constrain the parame- 
ters even further. Our conclusions are given in the final section. 

2. Observations and Data Reduction 

With the 10 m Keck I Telescope at the Mauna Kea Observatory 
two hundred high-resolution spectra were obtained by N. 
Przybilla in half a night in July 2004, using the High Resolution 
Echelle Spectrometer (HIRES; Vogt et al. [D^ . The spectra 
covered a wavelength range of 4 200 A - 6 600 A with a few 
small gaps at a resolution of 0.1 A and exposure times of 20 s 
each. The data were reduced using the ESO-MIDAS package. 
Bias and flatfield corrections were applied. The echelle orders 
were extracted and sky background was subtracted. The wave- 
length calibration was done separately for every single order 
All spectra were corrected to the heliocentric frame of refer- 
ence. To perform model atmosphere fits all spectra were radial 
velocity corrected to rest wavelength and a coadded spectrum 
was generated separately for every echelle order 

Two spectra were taken with the ESO Very Large 
Telescope UT2 (Kueyen) equipped with the UV- Visual Echelle 
Spectrograph (UVES) in spring 2001 (Dekker et al. 2000 1. The 
spectra cover a wavelength range of 3300 - 6650 A with two 
small gaps at a resolution of 0.1 A. The exposure time was 300 s 
each. The spectra were reduced with a pipeline developed by C. 
Karl using the ESO-MIDAS package and based on the imple- 



mented UVES reduction pipeline (Karl l20(J?l . Additional ob- 
servations were obtained with the 2.2 m Telescope at the Calar 
Alto Observatory in July 2004. The Calar Alto Faint Object 
Spectrograph (CAFOS) was used to obtain 150 spectra cover- 
ing a wavelength range of 3 600 A - 6 200 A with 5 A resolution 
and an exposure time of 180 s each. The data were reduced and 
coadded in analogy to the HIRES spectra. 

3. Radial Velocities and Spectroscopic Orbit 

Because the S/N ratio of the individual HIRES spectra was 
very low, only the Ha and H/3 lines could be used for deter- 
mining the radial velocity hy cross correlation with a model 
spectrum at rest wavelength. To improve the accuracy, the re- 
sulting radial velocity curve was combined with all available 



radial veloc ity dat a of KPD 1930+2752 (Maxted et al. 120001 
Woolf et al. 120021 Orosz 2000 priv. comm.) covering a time- 
span of four years. A sine curve was fitted to these 2 900 data 
points using a x' minimizing method and the power spectrum 
was generated (FITRV and FITPOW routines by H. Drechsel; 
Napiwotzki et al. 120011 . The sine curve fit is excellent for 
all different datasets and no period change could be detected 
over the whole timebase (Figs. ^ & |2j- The orbital parame- 
ters were measured with unprecedented accuracy and are con- 
sistent with prior measurements by Maxted et al. (.2000^ and 

1 kms"' 



Woolf et al. ( f2002b : r(Ha) = 5 + 1 km s ', = 341 + 1 km s"', 
P = 0.0950933 + 0.0000015 d. The residuals do not show any 
signs of eccentricity. 

No radial velocity variations due to pulsations could be 
found either. This result is consistent with the non-detection 
of RV variations by Woolf et al. (2002i. The high orbital RV- 
variation of the system enforces short exposure times in order 
to prevent orbital smearing effects. Considering the fact that our 
spectra were obtained with the largest optical telescope in the 
world, it can be doubted that such a detection would be feasible 
for this object with present-day instrumentation. 

4. Atmospheric models and synthetic spectra 

Up to now spectra of sdB stars were analysed either from metal 
line blanketed LTE model atmospheres or from NLTE model 
atmospheres neglecting metal line blanketing altogether As 
pointed out by Heber et al. 12000 1 and Heber & Edelmann 
(2004) systematic differences between these two approaches 
are present. Most importantly the gravity scale differs by about 
0.06 dex. O' Toole & Heber (i2006 ) found very high abundances 
of heavy elements in pulsating sdB stars of similar temperature 
as KPD 1930+2752 and showed that model atmospheres with 
supersolar metallicity (10 x solar) reproduced the optical spec- 
tra much better than solar metalhcity models did. Hence in the 
analysis of the spectra of KPD 1930+2752 we used their model 
grids for both solar and 10 times solar metallicity. 

As the gravity is of utmost importance for our analysis we 
also calculated new grids of models and synthetic spectra to 
account for NLTE effects and metal line blanketing simulta- 
neously. Since the temperature/density structure of an sdB at- 
mosphere is only slightly affected by NLTE effects (if at all), 
the LTE approximation is valid to this end. NLTE effects may 
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Fig. 2. Radial velocity curves for the years 2000 and 2001. Times are given in days. The corresponding HJD has to be added. 
Data are from Maxted et al. (2000: HJD = 2451651); Woolf et al. (2002 HJD = 2451716/2451717); Orosz (priv comm.: HJD = 
2451769/2451770), and from UVES spectra (HJD = 2452033/2452065). The soUd Hne is the fitted sine curve to all data points 
over a timespan of four years. The residuals in the lower parts show no signs of eccentricity. The error bars in the curve of Woolf 
et al. have been skipped, because they are very large due to the low S/N ratio of the spectra. 



become more important for the line formation of Balmer and 
helium lines. Therefore, we chose a "hybrid" approach by cal- 
culating the temperature/density stratification from metal line 
blanketed LTE model atmospheres and then performing line 
formation calculations for hydrogen and helium allowing for 
departures from LTE (Przybilla et al. 2006 1. We used both 
of the Kurucz codes, ATLAS9 and ATLAS 12, to calculate 
metal line blanketed model atmospheres (Kurucz .1993 1996j . 



The coupled statistical equilibrium and radiative transfer equa- 
tions were solved and spectrum synthesis with refined line- 
broadening theories was performed using Detae. and Surface 
(Giddings 1981. Butler & Giddings 1985). Both codes have 
undergone major revisions and improvements over the past few 
years. State-of-the-art NLTE model atoms for hydrogen and he- 
lium are utilised (Przybilla & Butler i2004i Przvbilla l2005l . 
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- KPD 1930+2752 

LTE (log m = 1) 




= 35183±85 K 
log g = 5,55±0,01 
log(n(He)/n(H)) = -1,47±0,010 




KPD 1930+2752 
LTE (log m = 0) 




= 35720 + 67 K 
log g = 5,67 + 0,01 
log(n(He)/n(H)) = -1,56 + 0,015 
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Fig. 3. Best fit with an LTE model of ten times solar metalHcity 
(model B of TableQ. 



Fig. 4. Fit with an LTE model grid of solar metalUcity (model 
AofTableU). 



The Stark broadening of hydrogen lines is important for the 
determination of temperature and gravity. Up to now all spec- 
tral analyses of sdB stars used synthetic Balmer line spectra 
based on the unified theory of Vidal, Cooper & Smith J1973I 
VCS) with improvements of Lemke (1997 ). Since new broad- 
ening tables for Balmer lines have become available (Stehle & 
Hutcheon I 19991 SH), we also calculated synthetic spectra from 
these tables to investigate systematic effects on the synthetic 
spectra. 

5. Atmospheric Parameter Determination 

The CAFOS spectra were preferred for the parameter deter- 
mination because they are of excellent S/N and the run of the 
continuum is smooth. The continuum of the Keck HIRES spec- 
trum suffers from a non-optimal rectification of the continuum, 
which rendered the analysis of the Balmer line wings difficult 
or unreliable. Therefore we restrained ourselves from deter- 
mining atmospheric parameters from the Keck HIRES spectra. 

The observed spectra were analysed with FITPROF by 
means of a fit (Napiwotzki 1999 ) ■ Rotational broadening 
was accounted for by choosing an appropriate value for the 
projected rotational velocity (see Section 6). 

Motivated by the discovery of strong log^ variations in 
the pulsating subdwarf star PG 1325-1-101 (Telting & 0stensen 
12004 *). this possibility was also taken into account. To search 
for changes of the stellar parameters due to pulsations or due 
to the ellipsoidal deformation of the subdwarf star, three to ten 
CAFOS spectra were coadded. The atmospheric parameters for 
every phase were determined. No changes could be detected. 
Then we combined all CAFOS spectra to form a very high S/N 
spectrum and performed a fit. The formal statistical errors were 
small, i.e. AT^s ^ 100K,Alogg * 0.005 dex, A log '-^ ^ 



0.005 dex. However, systematic errors due to inaccuracies in 
the model atmospheres and synthetic spectra may be much 
larger. To investigate this we used the different grids of models 
described in Section 0] to determine the atmospheric parame- 
ters. 

Table 1. Atmospheric Parameters of KPD 1930-H2752 derived 
from model atmosphere fits. A12=Kurucz ATLAS 12 code, 
VCS=Vidal, Cooper & Smith iW7^ . SH=Stehle & Hutcheon 
l fl999j . If^line formation. 





model 


[m/H] 




log,? 


n(H) 


A 


LTE VCS 





35 720 K 


5.67 


-1.56 


B 


LTE VCS 


+ 1 


35 183 K 


5.55 


-1.47 


C 


NLTE (If) VCS 


+ 1 


35 353 K 


5.61 


-1.50 


D 


NLTE (If) SH 


+ 1 


35 212 K 


5.67 


-1.51 


E 


NLTE (If) SH A12 


+ 1 


35 712 K 


5.67 


-1.58 




adopted 




35 200 K 


5.61 


-1.50 








±500 


±0.06 


±0.02 



We started the analysis by using the metal line blanketed 



LTE model atmosphere grids of O'Toole & Heber J2006> . 
Using solar metallicity models a simultaneous fit of the hy- 
drogen and some helium lines (Hei 4026A, Hen 5412A and 
in particular Hen 4686A) was not possible (Fit A in TableQ. 
This so called helium problem already occurred during other 
analyses of pulsating subdwarfs (Heber et al. 2000. Edelmann 
|2003l. The analysis of HST-UV spectra of three sdB stars 
with similar Tgff as KPD 1930-1-2752 revealed supersolar abun- 
dances of the iron group elements (O'Toole & Heber .2006j . 
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Those stars also displayed the optical He ionisation prob- 
lem. Using more appropriate metal -rich models (10 times so- 
lar metallicity) the problem could be remedied because of a 
modified atmospheric temperature structure due to significantly 
increased line blanketing (see also Heber et al. 2006 1. The 
abundances of the iron group elements could not been deter- 
mined for KPD 1930+2752 due to its high rotational velocity, 
which causes a strong broadening of all metal lines. Because 
of the similarity of its atmospheric parameters to those of the 
stars studied by O' Toole and Heber ( 2006 J . we adopted high- 
metallicity models as well and, indeed, the fit improved (Fit B 
inTablenpigs.m&Eli. 

In a second step we used the new grid of hybrid models 
that account for both metal line blanketing and departures from 
LTE using ATLAS9. We chose supersolar metallicity because 
of the experience with the LTE analysis described above. The 
quaUty of this fit (labelled C in Tabled is comparable to that 
from the LTE analysis. 

Then we checked the influence of Balmer line broadening 
by using synthetic spectra calculated from ATLAS9 models, 10 
times solar metallicity with the hybrid approach and SH broad- 
ening tables instead of VCS tables (Fit D in Table [Q. Finally 
ATLAS9 was replaced by ATLAS 12 and the analysis repeated 
(Fit E in Table using the same approach as in the previous 
step. 

The resulting parameters from the different grids are sum- 
marized in Tabled As can be seen, systematic effects for the ef- 
fective temperature are small, as models A to E differ by 530 K 
or less. The gravities, however, differ by as much as 0.12 dex. 
The "hybrid" models yield a gravity 0.06 dex higher than the 
LTE models for supersolar metal content. The replacement of 
the VCS tables by SH tables systematically increases gravity 
by another 0.06 dex. It is possible that part of the shift is ex- 
plained by differences in table organisation and interpolation 



("Lemke ri997t . Note in passing that the solar metallicity LTE 
model (VCS broadening tables), i.e. the one of the least so- 
phistication, yields almost identical T^ff and gravity as the most 
sophisticated model grid (supersolar ATLAS 12 plus NLTE line 
formation plus SH line broadening). However, the fit of the ob- 
served spectrum from the latter grid is superior. It is impossible 
to judge which of the models in Table^is preferable. The qual- 
ity of fit lends support for synthetic spectra calculated from su- 
persolar metallicity models. We therefore dismissed fit A (cal- 
culated for solar metallicity) and adopted the mean values of 
the other fits: = 35 200 ±500 K and log § = 5.61 +0.06 dex. 
For sake of completeness we mention that the helium abun- 
dance differs by 0. 11 dex. 

It is important to note that the effective temperature is sig- 
nificantly higher (about 2 000K) than previously derived in 
other investigations (Billeres et al. 2000) irrespective of the 
choice of model atmosphere. Thus, KPD 1930+2752 is situ- 
ated at the hot edge of the instability strip in the T^ff - log g 
diagram (Charpinet et al. ll996t . 




UVES 



\tff He II 4686 



HIRES 



A(A) 




A{A) 



Fig. 5. Best fits of Viot sin / to some lines of the HIRES and 
UVES spectra. The atmospheric parameters were fixed to the 
values derived from the CAFOS spectra. LTE models with 
NLTE line formation and ten times solar metallicity were used 
(model grid C, Table [TJ. 

KPD 1930+2752 as accurately as possible. For this purpose, 
a high S/N spectrum was constructed from the 200 HIRES 
spectra. Beforehand the individual spectra were shifted to rest 
wavelength. The median was calculated in order to filter cos- 
mic s. 

The projected rotational velocity was measured by convolv- 
ing a synthetic spectrum calculated from the best fit model 
atmosphere with a rotational broadening ellipse for appropri- 
ate v'lot sin /. The FITSB2 routine of R. Napiwotzki was used, 
which fits the projected rotational velocity and performes an er- 
ror calculation based on a bootstrapping algorithm (Napiwotzki 
et al. 2004 1. The lines Ha, Hy3, Hy, Heii 4686A, He i 4472 A, 
He 1 4922 A, and He 1 5016 A were used simultaneously for this 
measurement. The result is Vyot sin / = 93.8 + 2.2 km s"'. 

Fortunately one of the two UVES spectra was taken near 
the maximum of the radial velocity curve (see Fig. |2] lower 
right corner). That means that despite the relatively long expo- 
sure time orbital smearing is negligible and the spectrum can 
be used for the Viot sin /-measurement, too. The rotational ve- 
locity was obtained as described above from the lines Ho- - H8, 
He II 4686 A, He i 4026 A, He i 4472 A, He i 4922 A, and He i 
5016 A. The measured Vrotsin; - 90.7 + 2.1 kms"' is consis- 
tent with the HIRES result. To get a smaller error margin we 
calculated the average of the two independent measurements 
I'l-otsin/ = 92.3 + 1.5 kms"' (see Fig.|5}. 



6. Projected Rotational Velocity 

The primary aim of the high-resolution time-series spec- 
troscopy was to measure the projected rotational velocity of 



7. Mass and Inclination 

KPD 1930+2752 is obviously affected by the gravitional forces 
of the companion, demonstrated by its ellipsoidal deformation. 
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Fig. 6. Inclination versus sdB mass. The dashed lines corre- 
spond to + Icr error. The upper curve indicates the minimum 
incUnation for an eclipse. Since sin / cannot exceed unity, a 
minimum mass for the sdB follows. 

Since the period of the photometric variations is exactly half 
the period of the radial orbit and the two components of the 
binary are close together the rotation of the sdB star is very 
likely tidally locked to the orbit. Having determined the gravity 
and projected rotational velocity, we have three equations at 
hand that constrain the system, with the sdB mass MsdB being 
the only free parameter Besides the mass function 



/(MsdB,M,on,p) = 



these are 



(MsdB + M,omp)^ 2nG 



R 



2nR 



(1) 



(2) 



(3) 



With log g obtained from the model atmosphere analysis, 
the radius of the star R was calculated using the standard mass- 
radius relation (Eq.|3}. Together with the orbital period of the 
system P and the projected rotational velocity Viotsini = ^lot sin / 
the inclination of the system sin / was derived for different val- 
ues of the sdB mass. Because the rotation is tidally locked to 
the orbit, the rotational period of the sdB equals the orbital 
period of the system. Therefore the absolute value of the ro- 
tational velocity could be calculated. From the measured pro- 
jected rotational velocity the inclination of the system was then 
derived (Eq.|3. With the sdB mass as free parameter, the mea- 
sured radial velocity semiamplitude /TjdB and orbital period P 
the mass function was solved numerically (Eq.^ to derive the 



mass of the companion Mcomp and calculate the total mass of 
the binary. The fact that sin / cannot exceed unity gave a lower 
limit for the mass of the sdB M^dB > 0.45 Mq (see Fig.|6ll. The 
errors were calculated with Gaussian error propagation. The 
error budget is dominated by the error in log g, which has to 
be estimated from the model atmosphere fit. Period P, radial 
velocity semiamplitude /TsdB and projected rotational velocity 
Viot sin i are measured to a precision more than sufficient for 
our purpose. This means that the accuracy of the results can- 
not be improved before better model grids become available. 
The higher value of \ogg which was obtained by using new 
line broadening tables indicates a shift of the mass minimum 
and therefore of the total binary mass towards higher values. 
As can be seen in Fig. the total mass of the system exceeds 
the Chandrasekhar limit for almost all assumptions of MsdB- 
If the companion is a white dwarf, its single mass has to be 
lower than the Chandrasekhar limit. This implies an upper limit 
MsdB ^ 0.64 Mq and a possible total mass range of MsdB+WD - 
1.33 +0.08 -2.04 + 0.14 Mo (Fig.|7}. A more massive invisible 
companion such as a neutron star or a black hole cannot yet 
be ruled out completely. In this case the subdwarf mass would 
exceed the possible mass range for this kind of stars (Han et 
al. 120021 12003> . The inclination angle of the system (Fig. |6} 
is close to 90°. KPD 1930+2752 could be an eclipsing binary 
(cf. the very similar binary KPD 0422+5421, Orosz & Wade 
J1999t V The detection of eclipses would rule out any object 
smaller than a white dwarf. If eclipses were detected the sdB 
mass would be constrained to 0.45 M© to 0.54 M© (see Fig.|6}. 
On the other had, if eclipses could be ruled out from high qual- 
ity light curves, the total mass of the system must be larger than 
1 .55 Mq, well above the Chandrasekhar limit. 

8. Constraints by photometry 

The light curve obtained by Billeres et al. J2000t shows ellip- 
soidal variations. We used this information to further constrain 
the parameters of the KPD 1930+2752 system. We employed 
the light curve synthesis and solution code MORO which is 
based on the model by Wilson & Devinney J1971I I. The de- 
tails of the Bamberg implementation are given by Drechsel et 



al. J1995> . The software uses a modified Roche model for light 
curve synthesis. It is capable to realistically simulate the distor- 
tions of the stars induced by the presence of a companion. 

For comparison of the observations with the synthetic light 
curves of MORO, it was necessary to prepare the data accord- 
ingly. Since the program does not include the effect of Doppler 
boosting arising from the high orbital velocity, this effect had 
to be corrected for. It increases the total flux by a factor of 
(1 - v{t)lcf', where v is the radial velocity at time t and c the 
speed of light. The boosting is counteracted by the Doppler 
shift, which reduces the effect by a factor of (1 - v(t)/c)^ (see 
Maxted et al. ,2000) . A resulting factor of (1 - v(t)/c) was ap- 
plied to the total flux. Furthermore, the light curve shows a dip 
during the maximum at orbital phase 0.25. This effect is cur- 
rently unexplained and cannot be modelled by MORO, there- 
fore the corresponding data points have been left out for the 
light curve analysis. Because KPD 1930+2752 is a pulsating 
subdwarf with a very rich spectrum of modes, the short period 
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Fig. 7. Total mass of the binary system as a function of the sdB 
mass with associated error margins (dashed lines) for the total 
mass. The horizontal line marks the Chandrasekhar limit. The 
dotted vertical line marks the point where the companion mass 
equals this limit. For higher sdB masses, the companion can- 
not be a white dwarf, but has to be a heavier object such as a 
neutron star (NS) or a black hole (BH). 
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Fig. 8. Light curve data taken from Billeres et al. J2000I I shown 
together with examples of synthetic light curves. A good fit 
(A^sdB ~ 0.47 Mq; solid line) is compared to light curves for 
MsdB ~ 0.45 Mq and MsdB ~ 0.55 M©. For these two cases the 
other parameters have been calculated according to Eqns.^to 



pulsations have been filtered out by Billeres et al. (I2000> . As 
this procedure is not trivial at all, hidden pulsations may Umit 
the accuracy of the results. One the other hand, the amplitude 
of the ellipsoidal variations is very sensitive to the system pa- 
rameters as can be seen in Fig.|8] 



Fig. 9. Detail of Fig.|6lwith light curve data superimposed. The 
full drawn line is the relation derived from spectroscopy. The 
dashed lines mark the errors in sin /. Each dot denotes the lo- 
cation of a parameter set for which a light curve has been cal- 
culated. The background shading corresponds to the quality of 
the fit - darker shading implying better Note that good so- 
lutions cluster slightly above the eclipse limit, indicating a high 
orbital inchnation near 80°. 



At first trial, light curves with different component masses 
and orbital inclinations were synthesized. The corresponding 
changes in the light curve shapes were significant enough to en- 
courage further investigation. A coarse grid of light curves was 
then synthesized which covered the parameter space between 
A^sdB = 0.45 M0 and MsdB = 1 .0 Mq. For sin / and R^^^ a range 
of values was calculated for each MsdB according to Eqns. [J 
to |3] These ranges covered the parameter space between the 
respective error limits derived from spectroscopy. Other rele- 
vant parameters for light curve synthesis like orbital separation 
a and the companion mass were derived from Kepler's laws 
of orbital motion. The radius of the companion was calculated 
from the mass-radius-relation for white dwarfs, or set to an in- 
finitesimal value for neutron star masses (in that case ruling out 
eclipse effects in the light curves). A value was calculated 
for each grid point with respect to the obsered light curve. A 
coarse analysis revealed a marked tendency of good fits to clus- 
ter in the sin / regions of marginal eclipse. According to Eqns. 
|2]and|3]this behavior strongly favored sdB masses close to the 
canonical value (Fig.llO>. Lower values of sin / and higher sdB 
masses could be ruled out quickly, since they led to quick dete- 
rioration of fit quality. 

In order to probe the relevant portion of the parameter 
space in greater detail, a very fine grid of synthetic light curves 
was constructed, containing more than 76000 parameter com- 
binations. It covered only that area of parameter space where 
promising solutions had been identified previously. The range 
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of sdB masses extended from 0.45 M© to 0.60 Mq. Indeed, a 
significant overlap between the spectroscopically determined 
parameter range and a set of very good light curve fits could 
be found (see Fig. Most notable is the coincidence of the 
best fits with the region of eclipse (/ x 80°). Fit quality deterio- 
rates significantly for lower inclinations. Probably this behavior 
stems from the light curve minimum at orbital phase 0, which 
shows signs of an eclipse effect and is therefore fit marginally 
better by synthetic light curves displaying the same effect. On 
the other hand, this marginal dip could also be a remnant from 
the removal of the intrinsic stellar pulsations. A better light 
curve is needed to decide whether the binary is eclipsing or 
not. 

As noted previously, the light curve fits, favoring an eclipse 
and therefore a high inclination, constrain the possible sys- 
tem mass considerably. However the only absolute parameter 
which can be determined from the light curves alone is the 
orbital inclination, as the results of light curve analyses are 
well known to be degenerate with respect to the system scale. 
Indeed, the sdB mass is poorly constrained by the light curve 
fit (see Fig.|5) as are the orbital separation and the radii of the 
stars. However, if we combine the results of the photometric 
analysis with the mass-inclination relation derived from spec- 
troscopy (Fig.|6|l, the sdB mass is constrained to a very narrow 
range of MsdB - 0.45 - 0.52 Mq, corresponding to a total mass 
of 1 .36 - 1 .48 Mq (see Fig.[TOJ. On the other hand, sdB masses 
higher than MsdB ~ 0.52 Mq can be discarded since in that part 
of the parameter space the spectroscopic mass relation does not 
have an overlap with good light curve fits within the eiTor mar- 
gins. A neutron star companion can therefore be ruled out. 



9. Conclusion 

The combination of an extensive set of 350 high- and low- 
resolution spectra, 2 700 radial velocity data points from prior 
observations and a reanalysis of the published light curve made 
it possible to determine the parameters of KPD 1930+2752 
with unprecedented accuracy. The companion is a heavy white 
dwarf. KPD 1930+2752 is the first subdwarf whose mass could 
be constrained in this way. With 0.45-0.52 Mq it is remarkebly 
close to the canonical mass of 0.5 Mq. It is worthwhile to note 
that masses of other sdB V derived by asteroseismology are also 
close to the canonical value (Brassard et al. 1200 II Charpinet et 
al. '2005 i,b). Furthermore Orosz & Wade (19991 determine a 
sdB mass of 0.51 + 0.05 Mq for the eclipsing sdB+WD binary 
KPD 0422+5421. 

The total mass and the merging time of the binary indicate that 
it is an excellent double degenerate candidate for an SN la pro- 
genitor. KPD 1930+2752 is one of only two double degen- 
erate systems, which fulfill these requirements (see Fig. 
Napiwotzki et al. 2005 1. While it fits into the DD scenario, 
KPD 1930+2752 might also evolve into an SN la as a single 
degenerate binary system. Since the time for merging due to 
gravitational wave radiation fmergei- ~ 2 • 10'^ yr is of the same 
order as the EHB life time fEHB ~ 2 ■ 10'^ yr, Roche lobe over- 
flow could occur well before the sdB becomes a white dwarf 
due to the shrinkage of the orbit. Recent calculations by Han 
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Fig. 10. Detail of Fig. with light curve data superimposed. 
The full drawn line is the relation derived from spectroscopy. 
The dashed lines mark the errors in total mass. Each dot de- 
notes the location of a parameter set for which a light curve has 
been calculated. The background shading coiTesponds to the 
quality of the fit - darker shading implying better 

and Podsiadlowski (priv. comm.) indicate that the mass trans- 
fer would be stable. 

Follow-up observations should be undertaken to measure 
an improved light curve and search for clear signs of eclipses. A 
promising alternative option is provided for asteroseismology. 
KPD 1930+2752 is the first pulsating subdwarf with a mass 
estimate from orbital analysis. It could be used to calibrate as- 
teroseismological models (see e.g. Charpinet et al. 2005 1. This 
is a very difficult task because the frequency spectrum is very 
complex and modelling is rendered difficult by the ellipsoidal 
deformation and the rapid rotation of KPD 1930+2752. 
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